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ABSTRACT 

Using mid-infrared and radio selection criteria, we pre-select a sample of candidate 
high-redshift type-2 quasars in the Subaru XMM-Newton Deep Field (SXDF). To 
filter out star burst contaminants, we use a bayesian method to fit the spectral energy 
distributions (SEDs) between 24-^m and B-band, obtain photometric redshifts, and 
identify the best candidates for high-z type-2 quasars. This leaves us with f 2 Zphot ^ 
1.7 type-2 quasar candidates in an area ~ 0.8 deg^, of which only two have secure 
X-ray detections. The two detected sources have estimated column densities A^h '^2 
& 3 X 10^"^ m~^, i.e. heavily obscured but Compton-thin quasars. Given the large 
bolometric luminosities and redshifts of the undetected objects, the lack of X-ray 
detections suggests extreme absorbing columns A^h ^ 10^* m~^ are typical. We have 
found evidence for a population of "Compton-thick" high-redshift type-2 quasars, at 
least comparable to, and probably larger than the type-1 quasar population, although 
spectroscopic confirmation of their AGN nature is important. 
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1 INTRODUCTION 

High-redshift quasars with column densities A^h 10^* 
are so heavily absorbed that they are barely detectable 
in even the most sensiti ve hard X-ray surveys (see, e.g., 
iBrandt fc Hasingeij |2005| ). and are known as "Compton- 
thick". While in the local Universe, 40% of active galactic 
nucle i (AGN) are found to be Compton-thick (R isalit i et al.. 
1 19991 ). the fraction of high-redshift quasars with such ab- 
sorbing columns is currently unknown. 

An alternative to (hard) X-ray selection is in the mid- 
infrared, where the obscuration due to dust becomes small. 
AGN are powerful mid-infrared emitters due to dust sur- 
rounding the accretion disk (the torus) reprocessing the UV 
and X-ray photons from the cen tral engine. AGN in visible 
in X-rays were indeed found bv iDonlev et al.l l|2005l '). who 
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used a mid-infrared and radio excess criterion to to select a 
sample of AGN. The majority of these AGN are at z 1 
and are better described as Seyfert-2s by virtue of their low 
luminosities. The authors found that while the sources de- 
tected in the X-ray were unlikely to be Compton-thick, the 
AGN not detected in the X-ray (~20% of their sample) 
could be Compton-thick (iVn^ 10^** m"^). The obscured 
and gravitationally-lensed quasar IRAS FS C 10214-1-4724, at 
z ~ 2 .3 also appears to be Compton-thick ([Alexander et al] 

|2005| V 

iPolletta et al] (120061 ) use two different selection crite- 
ria. X-ray and infrared, to look for Compton-thick AGN in 
an area of 0.6 deg2. Their X -ray selection, with a flux limit 
~10"" W m"^ in the 2.5-8 keV band, finds 5 Compton- 
thick AGN. Of these 2 are spectroscopically confirmed high- 
redshift quasars. Of the infrared-selected sources, the strict 
SED criteria imposed on the sources guarantee AGN, but are 
likely to exclude heavily obscured [A-v <; 20) high-redshift 
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{z ^ 1) sources which do not show power-law mid-infrared 
spectral energy distributions (SEDs), but have significant 
contributions from stellar light at 3.6 and 4.5 ^m. It is not 
suprising, therefore, that they find a smaller Compton-thick 
fraction (~10%) than expected in their infrared-selected 
sample, as it is probably biased against heavily obscured 
AGN. At high luminosities and redshift, Alonso-Herrero et 
al. (2006, AH06) find ~50% of their sources are undetected 
in X-rays (although not necessarily Compton-thick). This 
sample is based on a power-law criterion, and could plausi- 
bly also be missing the most heavil y obscured sources. 

iMartmez-Sansigre et all l|2005h . hereafter MS05, found 
a population of high-redshift type-2 quasars at least as nu- 
merous as the unobscured (type-1) population and which 
possibly outnumber the type-Is by ~2-3:l. The most re- 
cent X-ray studies agree with a large obscured fraction but 
disagree on the r elative numb ers of Compton-thick sources: 
iDwellv fc P"ag3 (j2006. DPOel ) find a ~3:1 ratio amongst 
Compton-thin type-2 and type-1 quasars, while Gilli, Co- 
mastri & Hasinger (2006) infer a ratio of Compton-thin 
type-2 to type-1 quasars, although the total type-2 to type-1 
ratio can increase to ~2:1 when Compton-thick sources are 
included. 

In this letter, we present the X-ray properties of a sam- 
ple of obscured (type-2) high-redshift quasars, selected from 
mid-infrared and radio data in the Subaru XMM-Newton 
Deep Field (SXDF). 



2 SAMPLE SELECTION AND DATASET 

We proceed to select a similar sample to that of MS05, 
in the SXDF, where deep X-ray data are available. In 
this work we decrease the lower radio flux density crite- 
rion, to increase the number of candidates. This is im- 
portant as here we have a smaller area (~0.8 deg^ as op- 
posed to 3.8 deg'^ in the MS05 sample). The selection cri- 
teria used here are 5*24 jim > 300 ^Jy, Ss.e ,jm ^ 45 /ijy, 
and 100 piJy 5i.4ghz ^ 2 mjy, which yields 38 can- 
didates. For a detailed discussion of these c riteria see 
iMartfnez-Sansigre et all ()2006a. hereafter MSOeh . 

As a brief summary, the two mid-infrared criteria are 
able to target z 2 type-2 quasars with Ay ^ 5, but will 
also allow z ^ 1 ultra-luminous infrared galaxies (ULIRGs) 
in the sample. These ULIRGs will have radio lum inosities 
following the FIR-radio correlation (ICondonlll992l ). and so 
by choosing a high-enough radio flux, one can cut out all 
but the most extreme starburst contaminants. The MS05 
criteria were carefully chosen to avoid starburst galaxies, 
and lowering the radio criterion has the disadvantage that 
less extreme starbursts are allowed in the sample. 

The mid-infrared data were obtained by using the 
SWIRE DR2 (Surace et al. 2005), which covers the SXDF, 
and has a flux density limit of ^^250 /xJy at 24 /xm (50"), and 
10 /iJy for both 3.6 and 4.5 /im (lOcr and 5 a respectively). 
The 1 .4 GHz catalogue used is described in lSimpson et all 
l|2006ll and comes from a B-array (~5x4 arcsec^ beam) VLA 
survey with a peak flux density limit of 100 /xJy (~5-8 a). 

The spectroscopic completeness of the SXDF are poor 
at this stage, so we decided to undertake this preliminary 
study using photometric redshifts (see Section 3). The op- 
tical data were obtained from the Subaru B,V,R,i' and z' 



imaging (Furusawa et al., in prep.). The near-infrared data 
(J and K) w ere taken from the UKIDSS ultra-deep survey 
(UDS) DRl (|Warren et al.ll2007l '). The X-ray data were ob- 
tained with the XMM-Newton observatory, which covered 
the SXDF field in 7 pointings, with exposures of 50 ks each, 
except for the central pointing, which was observed for 100 
ks (Ueda et al., in prep.). 



3 PHOTOMETRIC REDSHIFTS AND 

FILTERING OUT ULIRG CONTAMINANTS 

In order to disentangle type-2 quasars from ULIRGs and ob- 
tain photometric redshifts, the SEDs of the candidates were 
fitted using a model consisting of three components (warm 
dust, galaxy light and blue light). The normalisations of the 
components were allowed to vary together with the redshift, 
z, although for a given Ay and z, the normalisation of the 
warm dust component was fixed to fit the 24 fim data. Two 
different models we re investigated, and the bayesian odds ra- 
tio (e.g. ISivia|[l996l ) was used to select between the models. 
These models differed only in their warm dust properties 
and consisted of a quasar (with dust) and a ULIRG. The 
redshift was given a flat prior between ^ z ^ 7. No real 
quasars have been detected at z ^ 7, justifying this cutoff. 

The stellar population was modelled by taking the SED 
of a, z — elliptical galaxy from Coleman, Wu & Weedman 
(1980). SEDs showing smaller 4000- A breaks were well fit- 
ted by the elliptical galaxy together with the blue compo- 
nent (see below). The elliptical was normalised to match the 
K-band luminosit y of an L* galaxy i n the local K-band lumi- 
nosity function of ICole et al.l (|200lh . The luminosity of the 
galaxy relative to the local LJ^ was then allowed to vary with 
a prior flat in log-space, between —1.3 ^logio(igai/LK)^ 
1.3. This choice brackets the reasonable range of host galax- 
ies for a powerful quasar accounting for passive evolution: 
from a very faint 0.05LJ; to a very bright 20LJ; galaxy. The 
choice of a prior flat in log-space means that the host galaxy 
is believed, a priori, to be as likely to lie between O.ILJ^ and 
LJ; as between LJ; and lOLJ^. This is more realistic than a 
prior flat in real space, which would imply a quasar is as 
likely to be hosted by an LJ; galaxy as by a lOLJ; one. 

The blue component, between 912 A and 5000 A, has 
a physical motivation as well as serving a practical pur- 
pose. It can represent the scattered light from the ob- 
scured quasar, blue light from young stars or a UV-upturn 
brighter than that of the template elliptical galaxy used (and 
therefore a smaller 4000- A break). The slope was therefore 
set to 5*1^ oc u~°'^, which is representative of both type- 
1 quasars and starforming galaxies. The normalisation of 
the blue component is given as a ratio of the luminosity 
to that of an galaxy at 5000 A, and was allowed to 
vary independently of the galaxy, with a flat prior between 
—2 $;logio(Lbiue/L* 2, spanning the range between 

a blue component fainter than an elliptical galaxy's UV- 
upturn, and as bright as the most powerful type-1 quasar. 
To avoid overfitting and losing the ability to discriminate 
between models, neither the stellar or the blue component 
are reddened by dust. 

The ULIRG dust component was m odelled using the 
models of ISiebenmorgen fc Kruegell (120071 ). hereafter SK07, 
only allowing variation of one parameter. Ay. For a given 
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(J2000) 
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/ 517 


ID052 


02 16 17.92 


-05 07 18.56 


1.90±0.05 


41. 2^ 


119 


0.85 


0.55 


13.0 


-10.0 


<0.7 


ID123 


02 19 28.76 


-05 09 08.81 


1.75±0.05 


40.2 


72 


0.65 


0.20 


5.2 


-9.7 


<2.92 


ID135 


02 19 01.89 


-05 11 14.22 




40.2 


2.2^ 


1.20 


0.15 


26.0 


-9.4 


<0.3 


ID142 


02 17 23.82 


-04 35 13.72 


4.05±0.05 


40.4 


6.7 


1.25 


0.45 


10.0 


-8.9 


<0.4 


ID147 


02 19 10.31 


-05 16 03.00 


1.80±0.05 


40.1 


35.9 


0.75 


-0.05 


28.6 


-13.4 


<0.3 


ID200 


02 18 15.71 


-05 05 10.34 


, 7^-1-0.05 


40.1 


72 


0.75 


0.50 


8.4 


-9.1 


<0.4 


ID249 


02 19 13.74 


-04 56 04.27 


1.75±0.05 


40.0 


35.9 


0.50 


0.00 


25.2 


-9.2 


<0.6 


ID342 


02 17 05.35 


-05 09 24.61 


3.85±0.05 


40.0 


9.0 


1.20 


0.30 


4.6 


-15.25 


<6.82 


ID345 


02 16 29.56 


-05 03 10.65 


2 00+°-"** 


39.8 


35.9-1 


1.00 


0.10 


10.4 


-11.1 


<0.6 


ID347 


02 18 09.64 


-05 18 42.42 


1.75±0.05 


40.2 


72 


0.55 


-0.30 


5.8 


-8.0 


1.7±0.95 


ID386 


02 17 25.11 


-05 16 17.27 


1 90+°°^ 


39.4 


17.9 


0.50 


-0.75 


6.2 


-11.5 


<0.1 


ID401 


02 16 23.02 


-05 08 06.76 


^•^'-'-0.12 


40.0 


35.9 


0.70 


0.40 


26.7 


-8.9 


3.4±1.35 



Table 1. Be s t-fit p arameters and X-ray fluxes for the high-redshift type-2 quasars. The names are from the radio catalogue of 
ISimpsoii et al.l 1I2OO6I ). and we define 5*17 = 10 W m ^ for convenience. The errors in Zphot a^r^ estimated, given our model SEDs, 
from the full-width half-maximum values from the marginalised probability distribution functions (PDFs) for Zphot- For most objects, 
the PDF has at least one secondary peak, always adjacent to the primary one, with peak value between 0.5 and 0.05 of the primary peak. 
The median value of logio[Lgai/i]^], 0.75, corresponds to a 5.6Lj^ ( where L*, is the break in th e local K-band luminosity function) galaxy 
at z ~ 2, which asuming passive evolution and using the models of lBruzual fc Charloti | |2003|) . would become a present-day galaxy. 
The limits quoted for X-ray non-detections are the fiux returned by the task EMLDETECT at the position of the source, plus the error 
quoted for this flux. This is approximately the same as a 2(t limit. ^This high value for Lbol is probably due to an overestimation of 
the Av. ^The reason for these high values for the fluxes is that the background noise is significantly higher in these two sources. ^This 
value of Av suggests a reddened type-1 quasar rather than a genuine type-2. *ID345 is a point source at K-band, which suggests the 
value of Ay is an overestimate. ^ID347 and ID401 have safe X-ray detections. The best-fitting SEDs for all 12 objects can be found at: 
[http:/ /www. mpia.de /homes / martinez/sxdftype2seds.pdf 



Ay and z, the bolometric luminosity of the ULIRG, Lboi was 
chosen to make the SED go as close as possible to the 24- 
jum data point. The values of Lboi are only restricted by the 
range available in the SK07 library: 10.1 ^ logio(-£/boi/i'0) 
^ 12.7 for our choice of parameters. Of the other parameters, 
the nuclear radius was fixed to 1 kpc, the ratio of luminosity 
of OB stars with hot spots to Ltot was fixed to 0.6, and 
the hydrogen number density was fixed to 10^° m'^ The 
models have a discreet set of values of the extinction (Ay = 
2.2, 4.5, 6.7,9, 17.9, 35.9, 72 and 119), so the prior for the 
Ay consists of a set of 5-functions (with equal probability) 
at these values. 

The mid-infrared SED for quasars was modelled using 
the lElvis et"aLl (|l994l '). hereafter E94, type-1 SE P, an d ob- 
scuring it with dust from the models of iPeil (|1992| ). For 
consistency with SK07, only Milky Way (MW)-type dust 
is used. This SED and dust model would allow us, in prin- 
ciple, to vary the values of Ay continuously, but to make 
the fitting procedure as fair as possible, we restricted the 
values of Ay to those present in the models of SK07 and 
assigned the same prior to them. This flat prior reflects our 
ignorance about the range in Ay, particularly for quasars. 
For each object, the quasar bolometric luminosity Lboi was 
not allowed to vary freely, but was fixed for a given Ay and 
z by the observed flux density at 24 fim. There is, however, 
no upper limit set on Lboi- 

For a given model, the likelihood of a given combina- 
tion of parameters is given by the n-dimensional probability 
density function 

p(data|model A, z, Ay, Lgai, ibiuo) = 



(1) 

(27r)"/2]Ja, 

i 

where 

n 

xi = ^ (model A, - data,) ^ /erf , (2) 

i 

model Ai is the flux predicted by model A, given the pa- 
rameters, over a given waveband, data; is the observed flux 
at that band, ai is the measurement error in that band and 
n is the number of bands. In the cases where ai was smaller 
than 10% of the flux density of the object, it was set to 10%. 
To treat non-detections at a particular band, we followed 
the following method: when the observed galaxy flux den- 
sity and the model fell below the flux density limit {5atand, 
where atand is the rms noise in that band), the band made no 
contribution to ■ When the model lay above the limit, the 
object was assigned a flux density Si and an error ai both 
equal to half of the flux density limit (so Si —ai — 2.5aband)- 
An undetected source would therefore only be lai away from 
zero flux density, and from the limit. This prescription was 
appropriate for our sources, selected to be faint and therefore 
close to the limits in most of the bands, to avoid excessively 

penalising nondetections in one band. 

To select between models, we follow [sivial (|l996l ) in cal- 
culating the odds ratio (OR), 

- p(d|B) p(B) ' 

where, for brevity, we refer to model A as A, model B as 
B, and data as d. We assume that models A and B are 
equally likely a priori, p{A)/p{B) = 1, implying ignorance 
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in the fraction of quasars c.f. starbursting galaxies. Thus the 
evidence for model y (where y £ {A, B}) is proportional to 
p(d|y), which is simply the likelihood (from eq 1) integrated 
over the parameter space spanned by the priors, i.e.: 

P(d|y) = 

J e"^y''^p(2:)p(Av)p(I/gai)p(ibiuo)d2;d^vdigaidI,biuc, (4) 

Xy is as defined in Equation 2, p(x) is the normalised prior 
probability distribution of parameter x and the integral is 
over the entire prior space. To make sure that our X-ray 
analysis did not include any ULIRGs, only objects with 
OR(quasar/ulirg)^ 100 were considered as quasars. This is 
equivalent to stating we are at least 99% confident of the 
quasar model over the ULIRG model. Note this criterion 
is ve ry similar to the so-called Jeffreys' criterion (jjeffrevd 
Il96lh . In addition, to make sure the fits were acceptable, 
only sources with evidence p(data|quasar)/(S'i7)^^ 10~^^ 
were accepted, where S17 = 10~^^ W m"'^ {Sn is defined 
to avoid numbers of order lO"'^"'' or smaller). Our range of 
evidences in Table 1, 1.3 x 10~^ ^p(data|quasar)/(S'i7)^^5; 
5.5 X 10~^^, corresponds approximately to values of reduced 
= 1.9 — 6.1. Finally, only sources with Zphot ^ 1-70 were 
kept, since Ly a is visible in optical spectroscopy, so the 
redshifts can be checked observationally, as long as this line 
is not obscured by dust on large scales (see MS06). Pre- 
liminary spectroscopy of the sample suggests our values of 
Zphot are in good agreement with spectroscopic redshifts (A. 
Martmez-Sansigre et al. in prep.). 

Of the 38 candidates in the sample, 12 (32%) follow 
our above criteria and are considered our best candidates 
for high-redshift type-2 quasars. Their best-fit parameters 
are quoted in Table [T] and we proceed to analyse the X-ray 
properties of these sources. Visual inspection of the optical 
images showed several cases where the blue light was unre- 
solved. This fits in well with scattered light from the quasar 
reaching us and justifies our choice of a "blue" component. 
Typically the best fitting normalisation for the blue compo- 
nent is ~5% of the normalisation of the quasar component, 
which is a reasonable scattering fraction. The K-band im- 
ages showed resolved sources in 11 cases, while ID345 is a 
point-source (suggesting a real Ay 5; a reddened type-1, 
see MS06). The SED and Ay of ID 135 (see supplementary 
Figures) also suggest a reddened quasar, rather than a gen- 
uine type-2, although since it is not point-like in K-band, the 
host galaxy probably contributes significantly to the near- 
infrared flux. 



4 X-RAY PROPERTIES OF THE TYPE-2 
QUASARS 

The type-2 quasars were cross-matched with the X-ray cata- 
logue of Ueda et al. (in prep.) with a flux limit of ~3x 10~^* 
W m~2 in the 2-12 keV band, but only 2 out of 12 (17%) 
were detected: ID347 and ID401. These objects have like- 
hhoods (EP_DET_ML) of 81 and 80 respectively, and are 
therefore clear detections (see DP06 for an explanation of 
the likelihood). 

To obtain meaningful limits for the undetected sources, 
these had their counts measured directly from the X-ray 



J_Ql3 ^Ql4 '^qVM"^ ^QlB 




^q13 ^qU ^q15 ^q18 ^q17 ^q18 
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Figure 1. Best fit spectral energy distribution (SED) of ID401. 
The top X-axis and right y-axis show the rest-frame frequencies 
and luminosities (in units of Lq) assuming the best-fit ^phot = 
1.90. The dashed black lines show (from left to right) the quasar 
(with extinction), galaxy and blue templates normalised by the 
values (from Table[T]l alongside the best fit to the data (red dots), 
the solid purple line shows the combined best-fit templates. The 
black dots and limits are the X-ray data in the 0.3-0.5, 0.5-1, 
1-2, 2-4.5 and 4.5-12 keV bands, from all three EPIC cameras. 
The coloured lines are the WF99 absorbed models, with values 
of logio(ArH/ m-2) = 25.00 (red), 27.25 (yellow), 27.50 (green), 
27.75 (cyan) and 28.50 (dark blue), with normalisation so the 
unobscured curve (25.00) goes through the highest energy X-ray 
point. The X-ray datapoints are best described by the curve with 
27.25 (yellow). For comparison, the E94 type-1 SED has been 
overplotted (dash-dotted line), with the same values of Lbol ^ud 
^phot as ID401 (from Table \T\ . The low value of the unobscured 
WF99 curve, compared to the E94 curve, shows this source has 
an intrinsically lower X-ray to mid-infrared ratio by a factor ~ 3, 
than the expectation of E94. 



image, using the XMM-Newton Science Analysis Software 
(SAS). The positions of the sources not detected in the 
XMM-Newton observations were added manually to the 
EBOXDETECT source lists of each observation. We then 
ran the SAS task EMLDETECT keeping the positions of the 
sources fixed. EMLDETECT performs a maximum likeli- 
hood fit on the distribution of observed counts of the sources 
previously extracted by the task EBOXDETECT. The fit 
uses the five diferent energy bands (in keV: 0.2-0.5,0.5-1,1- 
2,2-4.5,4.5-12) and the three EPIC cameras (Ml,M2,pn) si- 
multaneously. This procedure allowed to fit the sources as 
multiple components with separate point-spread functions 
(PSFs), and sources with no evidence of any X-ray detec- 
tion were still fitted. 

This procedure allowed us to obtain fits to the counts for 
all sources in the 5 EPIC bands, and these were converted to 
fluxes in each band, assuming a photon index F = 1.7 and an 
obscuring column of A''h = 3 x 10^^ m~^ (the default SAS 
values). These are not the most likely values for obscured 
quasars, so we estimate the effects of using more representa- 
tive values. Assuming instead F — 1.9, A^h = 3 x lO^'^ m~^ 
at 2: = 2 we estimate the change in fluxes to be only -1-9%, 
+3% and -2% in the 1-2, 2-4.5 and 4.5-12 keV bands, and 
therefore negligible compared to the large uncertainties in 
the count rates. All the 10 sources undetected in the Ueda 
et al. catalogue were found to have fitted values consistent 
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with background noise, even in the hard band (see Table 1). 
We therefore find that 83% of our sources classified as high- 
redshift type-2 quasars are undetected in the X-ray image, 
even down to a flux limit of ~3xl0~^* W m~^ in the 2-12 
keV band (the exact limit varies across the X-ray image). 



5 DISCUSSION 

We have found that 10 out of 12 of our sources classified 
as type-2 quasars, are undetected in an X-ray image with 
an approximate flux limit ~3xl0~^* W m~^ (in the 2-12 
keV band). For the two detections, we u se the photometric 
redshifts and the Monte-Carlo models of IWilman fc FabianI 
||1999. WF99I ) to estimate the absorbing columns. The esti- 
mates for ID347 and ID401 are logio(iVH/m^) = 27.50 and 
27.25 respectively (although ID347 requires a black body at 
10^ K in addition to the WF99 spectrum, to fit an upturn in 
the soft X-rays). For both these objects, the intrinsic X-ray 
to mid-infrared ratio seems to be slightly lower than the E94 
expectation by a factor ~2-5. 

For the undetected objects, it is not possible to esti- 
mate values of A''h. However, as Figure 1 shows, the sensi- 
tivity of the XMM-Newton observations and the "negative 
K-correction" in X-rays mean that our X-ray data are sensi- 
tive to heavily obscured Compton-thin quasars, even if, like 
ID347 and ID401, at X-ray energies they are intrinsically 
weaker than our first expectation from the median SED of 
E94. Only when logio(iVH/m^) ^ 28.50, do sources with 
^boi ~ 10^° W become too faint to be detected in the en- 
ergy bands covering 2-12 keV (rest-frame ~6-36 keV). 

To get a handle of the possible values of A^h, in Fig- 
ure [2] we show fiducial tracks for a model quasar (see Figure 
caption for details). From this figure, the two detected ob- 
jects in X-rays and two of the non-detections have values 
consistent with being heavily absorbed but Compton-thin, 
while the other 8 objects undetected in X-rays are likely to 
be Compton-thick. Spectroscopic confirmation of the AGN 
nature of our sources is now required. 

Inspection of Figure 2 shows most of our objects pop- 
ulate a region of the S'2-i2kcV versus S2A iim plane that is 
lacking in sources in previous work (e.g. AH06). The rea- 
son for this is that our selection criteria focus in on objects 
covering a fairly narrow range in redshift and bolometric 
luminosity over a fairly large sky area. Most previous stud- 
ies (e.g. AH[06) cover much larger ranges in z and luminosity 
with few objects like ours expected because of small sky area 
coverage and mid-IR selection techniques (such as mid-IR- 
power-law selection) that bias against the most obscured 
objects. 

Following MS05, we model the expected number of 
type-1 quasars following our 24- /im and radio flux density 
criteria, at z ^ 1.70, and in 0.8 deg^. We predict 2±0.8 type- 
1 quasars following these criteria, while we have 12 type-2 
quasars of which 8 are likely Compton-thick. Assuming the 
error in the modelled number of type-1 quasars, and Pois- 
son errors for the numbers of type-2 quasars, the type-2 to 
type-1 ratio (in the 68% confidence interval) is in the range 
3.0-12.9 while the Compton-thick-type-2 to type-1 ratio is 
in the range 1.8-9.0. 

Although we are clearly suffering from problems due 
to small-number-statistics, the implication of our study is 



10"" 



Vf 



A^= 35.9 
log,„«„=27.50 

logioA'H=29.g5 



lO-* 



10* 



Figure 2. Hard X-ray flux (52— I2kev) versus 24-iim flux density 
(■S'24 fim) as a function of redshift for a model obscured quasar. 
This quasar has been chosen to have logio(ijbol/W^) = 40.1 and 
Ay = 35.9, which are both the median values from Table 1. The 
quasar has th e in t rinsic SED of E94, with obscuring dust from 
the models of iPeil lll992h . We use the absorbed X-ray spectra of 
WF99, with four different absorbing columns, logio(AfH/m~^) = 
26.75, 27.50, 28.25 and 29.25. This corresponds to a range in 
gas-to-dust ratio of A^hMv = 1-9 X 10^^ to 5 X 10^^ m'^ (the 
lower value being that of the Milky Way, the upper value being 
much larger but reasonable, see MS06 and Watanabe et al., 2004). 
The green and blue tracks are obscured but Compton-thin, the 
red track is mildly Compton-thick, while the yellow line is heavily 
Compton-thick. The black dots mark the values at particular red- 
shifts. Changing the Av to any other value of Table 1 would move 
the tracks horizontally only. The 12 type-2 quasars in our sample 
are overplotted in blue. The vertical dotted line is the 300 /ijy 
flux density limit imposed on our sample, the horizontal dotted 
line is the approximate flux limit of the X-ray image, ~ 3 X 10~^* 
W m~^. For our two X-ray detections, the independent estimates 
of logio ( Afjj /m'^ ) are ~ 27.50 so at flrst one expects them to fol- 
low the blue track. That these two objects lie signiflcantly below 
the blue track is not a cause for serious worry as it assumes a 
fiducial unobscured QSO SED that is subject to population vari- 
ance. To make the objects ID347 and ID401 consistent with the 
track would mean they have intrinisc X-ray to mid-infrared ra- 
tios ~ 3x lower than those of the assumed SED (E94, see also 
Figure 1). This ratio is consistent with the ~ la spread of the 
population used in E94. While the two non-detections with high 
limits are approximately consistent with ID347 and ID401 and 
the Compton-thin tracks, the remaining 8 sources are well below 
the expectations, suggesting they are Compton-thick. 



that, at high redshift, Compton-thick quasars may be the 
dominant sub-population of quasars. Although such objects 
would be essentially absent from even moderately hard X- 
ray surveys (energy ~2-10 keV), and their contribution to 
the hard X-ray background (energy J5IO keV) is diluted by 
their large distances, they could clearly represent a vital part 
of the accretion history of black holes. 

Our work to date (MS05, MS06) has included radio as 
well as mid-IR selection criteria. There is therefore a residual 
worry that the properties of the objects we have studied are 
infl uenced in some way by t h e prese nce of weak radio jets, 
and iMartfnez-Sansigre et al.l l|2006bl ) have shown that their 
radio properties are consistent with these objects having de- 
veloped large-scale FRI-like jets. In principle, only mid-IR 
selection criteria are needed to find high-redshift quasars, al- 
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though the difficulty Ues in filtering out the starbursts from 
the quasars. 
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